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Abstract 

One of the key elements of Feynman’s formulation of non-relativistic quantum mechanics is a 
so-called Feynman path integral. It plays an important role in the theory, but it appears as a 
postulate based on intuition rather than a well-defined object. All previous attempts to supply 
Feynman’s theory with rigorous mathematics have not been satisfactory. The difficulty comes from 
a need to define a measure on the infinite dimensional space of paths and to create an integral that 
would possess all of the properties requested by Feynman. 

In the present paper, we consider a new approach to defining the Feynman’s path integral, based 
on the theory developed by P. Muldowney. Muldowney uses the Henstock integration technique, 
and non-absolute integrability of the Fresnel integrals in order to obtain a representation of the 
Feynman’s path integral as a functional. This approach offers a mathematically rigorous definition 
supporting Feynman’s intuitive derivations. But in his work, Muldowney gives only local in space- 
time solutions. A physical solution to the non-relativistic Schrodinger equation must be global, 
and it must be given in the form of a unitary one-parameter group in L 2 (M n ). The purpose of 
this paper is to show that one-dimensional Muldowney’s local solutions may be extended to yield 
a global solution. Moreover, the global extension can be represented by a unitary one-parameter 
group acting in L 2 (M). 

PACS numbers: 02.30.Cj, 02.30.Sa, 02.30.Tb, 03.65.-w. 

Keywords: Schrodinger equation, Feynman’s path integral, Henstock integral, path-space average, global 
solution, unitary one-parameter group. 
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1. INTRODUCTION 

1.1 Motivation 

We consider three problems: (I) the definition of the Feynman path integral, (II) an 
approach to Feynman’s question via the Henstock integral, introduced by Muldowney jl|; 
and (III) a representation of the solution to the non-relativistic Schrodinger equation as a 
path-space average. Our emphasis here is on the word “average,” and we make this precise 
below. 
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In more detail, our first question (I) entails a closer look at Feynman’s suggestion for a 
representation of the solution to the non-relativistic Schrodinger equation as a path-space 
integral { 2 - 4 1. By the Schrodinger equation, we mean 

di’ _ 1 u , 

dt ih H 

^(0, x) = ^ 0 {x), 

where i/jq G F 2 (M n ), with n- fixed; and the Schrodinger operator has the form 

V(x); 


( 1 ) 


h 2 

H = - A x 

2m 


or setting the constant h — 1, we have 

dil> 


dt 


= 1 


— A x — V (x) 
2 m 




( 2 ) 


(3) 


^(0, x) = i/> 0 (x). 

By a Feynman integral, we mean a formal representation 

^(x, t) — ’’const” [ 


(4) 


where is the space of all paths starting at x, to : R 9 t (-)■ c u(t) G M n , cj(0) = x. 

In the more traditional approach in earlier literature to this problem (see e.g., [5|), the 
effort by many authors has concentrated on answering Feynman’s question with variations of 
what goes by the name “Feynman integral”, and hence a Feynman measure. It is well known, 
that a mathematical rigorous version of Feynman measure is fraught with difficulties, if not 
contradictions. While the relevant literature following this idea, and aiming for a Feynman 
measure, is vast, we mention here only that it is known that there is no sigma-additive 
positive candidate of Feynman measure; that is “the Feynman measure” is not positive 
measure on path-space. Also, the action term S(oj,t) in (J4J) is given by the time integral 
S(u>,t) = fg (y^ 2 ~ ^( a; (' s ))) ds, and is interpreted only as a formal limit (n —» 00 ) of 

2 


*SV) (xqj 


Nelson 


, in fl, 



V^,) 1 * 


(5) 


suggested the following approach to derive the Feyman path integral (JU). Fix 

1 \ 1/2 


m and V, and set operators: 

K^oix) = 


2n it J 

M t 'i/j 0 (x) = e~^ tv( - x ^o(x 


e2 * |x yl ipo(y)dy, 


R 


( 6 ) 

(7) 
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Then if H in ([2]) is essentially selfadjoint, Trotter’s formula converges in the strong operator 
topology on £>(L 2 (M n )) to the solution (J3J) of the Schrodinger equation dTJ) : 




U^o 

exp I it 



V’o 


lirn (K t/n M t/n ) n ^ 0 . 

n—too v / 


( 8 ) 


1.2 Our approach 


Our approach (II) (see theorem [H] sec 4.3) serves to reconcile two physical principles. The 
first one (i) is that time in diffusion is not reversible (the Second Law of Thermodynam¬ 
ics). Mathematically, the corresponding Feynman-Kac formula is a diffusion equation. The 
second principle (ii) is at the foundation of the meaning of the Schrodinger equation. The 
Schrodinger equation is an equation for the dynamics of wave-functions (representing quan¬ 
tum states). Energy is a conserved quantity, and the Schrodinger equation is time reversible; 
hence represented by a unitary one-parameter group acting in L 2 (M n ). Our approach is mo¬ 
tivated by a diffusion equation; but, to reconcile the two physical principles, (i) and (ii), we 
must introduce complex “transition-probabilities.” 

Our suggestion here is to give a representation of the solution to the Schrodinger equation, 
not as a path-space integral, but rather as a linear functional on functions defined on path- 
space. This approach will still yield the space-time solution to the Schrodinger equation, 
and as an average over path-space (3), but just not with respect to a positive sigma-additive 
measure. Still our approach has attractive features which are required of a solution to 
Feynman’s question; perhaps in the form Feynman should have asked it. Our solution is 
based on the Henstock integral, introduced by Muldowney 1]. 

When the machinery of the Henstock integral is applied to the problem, we do get a 
solution to the Schrodinger equation. However Muldowney’s work offers only local solution, 
so a solution defined only locally in space-time. Indeed this solution is purely local. But of 
course a physical solution to the non-relativistic Schrodinger equation must be global, and 
it must be given in the form of a unitary one-parameter group in L 2 (R”). Our purpose here 
is to show that a field of local solutions may be merged, and extended (in fact, uniquely) 
to yield a unitary one-parameter group U tl where t stands for time, acting in L 2 (W n ) and 
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producing now a global solution to the Schrodinger equation, with energy-conservation, and 
with time-reversal. 


Now there are alternative approaches: some based on analytic continuation in time, and 
others in mass, thereby making a connection to Wiener measure; and the Feynman-Kac 
formula for the solution of the inhomogeneous heat equation. Since Wiener measure is 
indeed well defined as a positive and sigma-additive measure on path-space, this approach 
has shown some promise, but it has also run up against some difficulties, see e.g., [5]. For 
example, even in favorable cases, there are analytic continuation from Wiener measure, but 
valid only for some values of mass parameter in the corresponding Schrodinger equation. 
Even so, the analytic continuation lacks the kind of direct physical interpretation implicit 
in Feynman’s original question. 


1.3 The Trotter Formula 


As always, the Hamiltonian operator for the Schrodinger equation is a non-commutative 
sum of a kinetic term, and a potential term. In general such a sum is only formally selfadjoint, 
but not selfadjoint in the sense of having a unique spectral resolution; the sum may have J. 
von Neumann defect indices different from (0, 0). Therefore, as a premise in our theorem, 
we must assume that this sum is indeed essentially selfadjoint as an operator in L 2 (M n ) (if 
there are N particles, then n = 3IV). Fortunately, there are already powerful results available 
(see e.g., 0 ) with realistic conditions on the potential which imply essential selfadjointness. 
And hence they apply to all physical potentials. With the tools mentioned above, we then 
show that the following two solutions agree: a local one and one based on a unitary one- 
parameter group. The Erst one, the local solution is based on an application of the Henstock 
integral, and the other, the global solution, is based on an approximation; in detail, a global 
time-space solution to the Schrodinger equation as a unitary one-parameter group U t , is 
achieved with the use of a Trotter-product formula approximation; so an infinite product 
representation; and with a Lagrangian representation in phase space. The convergence of 
the Trotter-product formula in turn requires that essentially selfadjointness holds for the 
Hamiltonian operator under consideration. 

The study of the Schrodinger operator and their realization as generators of one-parameter 
groups of unitary operators has been studied over decades. While these approaches vary, they 
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arc all different from ours. The questions addressed differ from ones of these inve stig ations. 


To give a sample, we mention the following papers, and the works cited there 
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23 ]. 


2 PRELIMINARIES 

2.1 Henstock Integral in Finite-Dimensional Space 

In this section, we will define the Henstock integral over a domain of a finite-dimensional 
space, R n . The set-up and most of the notation is taken from ]lj. 

Definition 1. A cell or cylindrical interval I ofW 1 is defined to be 

I = I(N) =fix-xl„, 

where N = {l,...,n} is the list of the dimensions; Ij, 1 < j < n, are one-dimensional 
intervals, that can be one of the following intervals: (—oo,a], (a, b\, (b, oo), (— 00 , 00 ). The 
collection of all cells in R n is denoted by I(R n ) = {/(IV)}. The union of a finite number of 
cells is called a figure and denoted be E. The collection of all figures in R n is denoted by 
E(R n ). 

To consider the Riemann sums for an integral over R n we associate with each cell /(IV) 
a sample point or a tag: x = (aq,... ,x n ) G R n . The sample point together with the cell 
form an associated pair ( x , /(IV)). We choose a sample point to be one of the vertices of the 
corresponding cell /(IV), that is for each 1 < j < n we have 

i) Xj = —00 if Ij = (— 00 , a], 

ii) Xj = a or Xj = b if Ij = (a, 6], 

iii) Xj = 00 if Ij = (6, 00 ), 

vi) Xj = —00 or Xj = 00 if Ij = (— 00 , 00 ). 

Since the tag points are allowed to take values of infinity, we need to introduce a real line 
with “points at infinity” adjoint. Let’s denote it by M = M U {— 00 , 00 }. 

Definition 2. A gauge in M n is a positive function 6 : M” —> (0, 00 ). An associated pair 
(. x,I(N )) is said to be 5-fine with respect to the given gauge 5(x), if for all 1 < j < n: 
(. Xj,Ij) is 5-fine in R, i.e. 

i) if I j = ( aj,bj\ then \Ij\ = bj — aj < 5(x), 
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ii) if Ij = (-00, a,] then a, < -jp^y, 
in) if Ij = (bj , 00) then bj > 

For the finite intervals Ij = bj\, the condition \I 3 \ < 8(x) means that in the limiting 
process the length of the interval Ij gets smaller as 8(x) becomes smaller. This condition 
is similar to the condition for choosing partitions in the definition of the proper Riemann 
integral on a bounded domain. The main difference is in the value for the bound. In the case 
of the Riemann integral, the bound is constant for each of the subintervals of a partition or 
the maximum of lengths of the subintervals is bounded by a constant. In the case of the 
Henstock integral, we replace the constant bound with a value of a gauge function S(x). The 
gauge function 5(x) can take different values depending on a particular sample point x and 
any particular properties of an integrand. This special feature of a gauge function makes 
some non-integrable functions in the Riemann or Lebesgue sense integrable in the sense of 
Henstock. 

In the case of unbounded intervals, Ij = (—00, a,] and Ij = (bj, 00), the corresponding 
conditions, a,j < — p^y and bj > pp, imply that as 5(— 00) and <5(oo) become smaller, 
the absolute value of an upper boundary for (— 00, af] and the value of a lower boundary 
for (bj, 00) become larger . Thus the unbounded intervals Ij = (—00,0,-] and Ij = (bj, 00) 
“shrink” in size. 

Definition 3. A partition in M n is a finite collection V of disjoint cells I(N ) whose union 
is M n . A division V in is a finite collection of associated point-cell pairs ( x , I(N )) whose 
cells I(N) form a partition ofW 1 . Given a gauge 5 in M n , a division V is called 5-fine if 
each of the pairs ( x , I(N )) EV is 5-fine. 

An integrand in M” is a function of associated points and cells. It is more intuitive to 
think about an integrand in the form of a product: f(x)h(I(N)). In this case a function 
h(I(N )) plays a role of a measure on the cells. But an integrand can be also given in a more 
general form h(x, I(N)). 

Definition 4. A function h(x, I(N )) is integrable on R n in the Henstock sense with integral 


a = 


h(x,I(N )), 


if, given e > 0 there exists a gauge 5 in M n so that, for each 5-fine division D$ ofMA, the 



corresponding Riemann sums satisfy 




a 


< £ 


or more explicitly 


\j2{Hx,I(N)):(x,I(N))eD s }-a 


< £. 


We defined the Henstock integral of a function with domain of integration being the 
whole space R n , but similarly the definition of the integral can be formulated for a figure E, 
a subset of R n . 

The definition of the Henstock integral relies on the fact that for every positive function 
5(x), a 6 - fine division T> of the domain of integration can be found. This result is due to 
the Belgium mathematician Pierre Cousin (1895) and is known in the literature as Cousin’s 
lemma or Cousin’s theorem. The proof in a one-dimensional case can be found, for example, 
24] or ]lj. In the finite-dimensional case of R n with n > 1, the result can be obtained by 


m 


applying the lemma in each dimension separately. 


2.2 Henstock Integral in Infinite-Dimensional Space 

In the previous section, we gave the definition of the Henstock integral in the finite¬ 
dimensional case. The gauge function S(x) and the 5-fine property of the point-cell pairs are 
essential for selecting appropriate partition in the definition of the Henstock integral in R n . 
For the case of infinite dimensions we will need to introduce more advanced kind of gauge 
and adjust gauge-fine property accordingly. In this section, we introduce necessary changes 
and then define the Henstock integral in the infinite-dimensional case. 

For the purpose of applying the Henstock theory to the Feynman path integral, the 
domain of integration that we will consider, is R T . R T is a product of T copies of R. By T 
we denote a real interval, a subset of (0, +oo). Usually, we assume T = (0, t] or T = (0, +oo). 
It would help our intuition if we would think of R r as of the space of all real-valued functions 
x(t) defined on T : 

R r := (x | x : T i->- R}. 

Note, R T includes both continuous and discontinuous functions, so it is not equal to the 
space of all paths, that are represented by continuous functions. 
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Let us denote by A f — Af(T) = {N = {t\,... ,t n } C T, n G N} a collection of all finite 

subsets of points in T. Suppose all the sets N = {fi,... ,t n } are ordered in the increasing 

manner: t x < t 2 < ■ ■ ■ < t n . 

Definition 5. A cell or cylindrical interval I ofR T is defined to be 

I = /[TV] = h x • • • x I n x R t \ n 

= h X • • • X I n X Y[ {K {i} | t e T\N} 

where N = {t\,... ,t n } G AfiT), Ij = I t for 1 < j < n. Similarly to the finite-dimensional 
case, Ij = I tj are the one-dimensional intervals of the form: (— 00 , a], (a, b\, ( 6 , 00 ) or 
(— 00 , 00 ). 

Note that to distinguish a finite-dimensional cell in R N from an infinite-dimensional one 
in M 7 , we use parenthesis in the notation in the finite-dimensional case: 

I(N ) = h x • • • x / n 

and brackets in the infinite-dimensional one. Thus a cell in M T can be written in the following 
form: 

I[N] = I(N) x M niV 

The collection of all cells in M 7 is denoted by I(M 7 ’) = {I[N} \ N e Af(T)}. A figure in M T 
is the union of a finite number of cells, denoted by E. The collection of all figures in M T is 
E(M t ). 

The form of functions that we will consider later, makes it appropriate to take a triple 
instead of a pair as a basic element of a partition. We consider triples with included set of 
restricted dimensions of the form (x, N, I[N]). The triple (x, N, I[N]) is called associated if 
I[N] — I\ x • • • x I n x R T \ N , N = {fx,..., t n } E Af(T) and x(N) = ( x(ti ),... ,x(t n )) is a 
vertex of I\ x • • • x I n , i.e. the corresponding finite-dimensional pair ((x(ti ) 7 ..., x(t n )), /(AT)) 
is an associated pair in R N . Note that if (x, N, I[N]) are associated and t G T\N , then 
x(t) G can be an arbitrary real number. 

Definition 6 . A gauge in 7 is a pair of mappings ( L,S ): 

L : t T -> Af(T) 
x (->• L(x) 


10 



and 


5 : R t x A f(T) -)■ (0, +oo) 

(x,N) S(x,N). 

With N = G MiT), an associated triple (x, N, I[N]) is 7 = ( L, 5)-fine or fine 

with respect to the gauge 7 = (L, 5 ) if 

1. L(x) C IV. 

2. For all tj G IV: (. x(tj),Ij) is 5-fine, i.e. 

i) if Ij = (a,j,bj\ then \Ij\ = bj — aj < 6 (x,N), 

H) = (-00,%] then aj < - S( ^ N) , 
in) if Ij = (bj,+ 00 ) then bj > g(+0 ^ N) . 

The first condition above means that the set L(x) is included in the dimensions where 
the cell I[N] is restricted. The second one requires that the finite dimensional point-cell pair 
((x(ti),..., x(t n )), I(N )) is fine with respect to gauge 5. 

Similarly to the finite-dimensional case, a division in M T is a finite collection D of point- 
ccll pairs (x,/[iV]) such that the corresponding triples (x, N, I[N]) are associated and the 
cells I[N] are disjoint with union Mfi. 

Definition 7. If a gauge 7 = (L, 5) is given, then a division D is 7 -fine, if each (x, N, /[iV])G 
D is 'y-fine. 

Definition 8. A function h of associated triples (x, N, I[N]) is integrable in in the 
Henstock sense with integral 

a — f h{x,N,I[N ]), 

if, given e > 0 there exists a gauge 7 so that, for each 'y-fine division D 1 of R T , the corre¬ 
sponding Riemann sums satisfy 

(D^J 2 h(x,N,I[N])-a <e 

or more explicitly 

I^ {h(x, N, J[iV]) : (x, N, I[N ]) G Dfi\ — a <e. 
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Note that in the case of R T , the rule 7 for selecting associated triples in Riemann sums has 
two conditions. We still use a positive function 5(x) to bound the lengths of the restricted 
edges of each cell I[N], but we also have an extra condition for the sets of restricted di¬ 
mensions N. In the finite-dimensional integral over R n , the set of restricted dimensions was 
always fixed. Now it becomes a variable. The condition requires, that the set of restricted 
dimensions N for each cell /[TV] in the partition, includes some minimal set of dimensions, 
given by a value of the gauge function L(x) at the corresponding sample point. Thus if we 
make 5{x) successively smaller and L{x) successively larger the cells in the corresponding 7 
- fine associated triples, will be “shrinking.” This is the idea in the limiting process of the 
Henstock integral in the infinite-dimensional space. 

Besides gauges, the specific feature of the definition of the Hensock integral is the im¬ 
plicit fact, that for a given gauge 7 , there exists a 7 - fine division Z ) 7 of R T , as in the 
finite-dimensional case. However, this time the proof is more technical and requires some 
extra tools from the theory of the Henstock integral. 


Theorem 1. Given any gauge 7 = (L,5) in R r , there exists a 7 - fine division D 1 of R r . 
Proof. See Theorem 4 in chapter 4 in [1], □ 


2.3 Elements of Henstock - Probability Theory 

Definition 9. A real- or complex-valued function F defined on the figures E(R T ) of R T is 
an additive cell function if F is finitely additive on disjoint figures. An additive cell function 
F is a probability distribution function if F( R T ) = 1. 

In classical probability theory a probability distribution function is also assumed to be 
non-negative. But for the case of the Feynman path integral, we drop this requirement 
and allow the distribution functions to be negative or complex-valued. This assumption 
makes our approach different from the classical one, when probability distributions are non¬ 
negative and correspond to the well-defined measures. Since we already know, that the 
classical approach does not guarantee satisfactory results for Feynman’s theory, our method 
based on including non-positive distribution functions, is justified. The motivation for our 
approach comes from a simple experiment with a beam of photons, and our argument lies 
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at the foundation of the wave-particle duality; i.e., the idea that elementary particles (such 
as photons-light quanta) behave simultaneously as particles and as waves. 

The experiment goes as follows: a laser beam is directed towards a particle detector- 
device. Using sensors we can measure a number of light particles that reach the detector 
with a certain probability. Since the particles also possess the wave properties, depending 
on their phase, the interaction between photons can result in cancellation, or intensification, 
of the beam. Thus a probability distribution function that measures the number of particles 
must then necessarily take on negative as well as positive values. This explains why, in 
quantum mechanics, it is reasonable to have a probability distribution function which attains 
complex values. 

Definition 10. If F is a probability distribution function onR T and if there exists a function 
f(x(N)) such that for each I G I(M T ) 

F(I) = Jj(x(N))\I{N]\ 

then f is a density function for F. Here by |J[iV]| = |/i x • • • x /„ x we understand 

the finite product 

\I[N]\ = \I 1 \-\I 2 \...\I n \. 


3 LOCAL SOLUTION TO THE SCHRODINGER EQUATION 


In his work Muldowney considers a general version of the Schrodinger equation in one- 
dimension of the form 


dfj(x,t) 1 d 2 fj(x,t) . . 

J iir 1+ Tc^F 1+cV ^ t ^ = 0 - 

Making the parameter c to be — | and the potential to depend only on space variable 
V(x,t) = V(x), we arrive to the Schrodinger equation in commonly used form with mass 
parameter taken to be 1 and h — 1: 

dfj(x, t ) 


dt 


= % 


121 _ l v(x) 

2 dx 2 2 V[ ' } 


ipfat). 


( 9 ) 


Below is the statement of the theorem as it appears in the work by Muldowney with the 
only change that the parameter c is taken to be — 
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Theorem 2. (Muldowney, 2012) Let and t > 0 be given. Suppose the function V(y) 

is continuous at y = x. If 

if{x,t) = E* t [U(x T -)\ 


exists in a neighborhood of (x,t) thenil)(x,t) satisfies the Schrodinger equation 

dtf(x, t ) 


dt 


= i 






( 10 ) 


Here E* t [U(xr-)\ is the marginal expectation given by 


E-JU(xt-)] = U(x t -)G(I[N -]). 


( 11 ) 


In this work, Xt represents a possible path of a particle; and is a real-valued function on 
time interval Xt '■ [0,£] —> R. By dehnition of the marginal expectation the paths that 
are considered in the integral (fill) are conditioned to have the right end-points fixed. The 
notation xt- is used to distinguish the path with a fixed right end from a path with a free 
right end. Thus the domain of integration for E* t \U(xT-)\ is given by 


R T = {x(t) | x : [0,i] —y R, x(0) = 0,a;(f) = x} 


and the integrands have the following definitions 

exp(—| f T _ V(x(t))dt), xt is continuous 



0 , otherwise. 


9 (x(JV-))|J(JV-)| 


g(x(N )) 


dx i... dx n -1 


#(0)=0,£(£)=:e 


—I 


2 ^ (tj tj- 1 ) 


i (^j 

e 2 


dxj. 


( 12 ) 


The function exp 


(Xj-Xj-1)* 


x(0)=0,x(t)=x 

under the integral sign is not Lebesgue integrable. How- 


^2 (tj tj— i) 

ever, we can use a closed contour in the complex plane to compute the improper Riemann 
integral for this function (see also [25-2^|). In one-dimension, the improper Riemann inte¬ 
gral is equivalent to the Henstock integral Q. Thus the Heustoch iu.egta. of the hVesnel 
integral with the simple expression in the exponent can be evaluated as follows 

f i x 2 , / 2 F 

/ e 2 dx i = \ — : . 

./ns V i 
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The method with a closed contour in the complex plane for one-dimensional case lies in the 
core of the proof for the result stated in the theorem below. Using the appropriate change 
of variables, the Riemann sum for the infinite dimensional integral can be computed as a 
product of one-dimensional integrals. We refer interested reader to [1] for details. 

Theorem 3. The Fresnel infinite-dimensional integrand is integrable in M T and 

[ G(I[N]) = 1. 

The function G (1121) is additive and has integral 1 over M r , thus it is the complex “tran¬ 
sition probability.” It represents the probability of a path to go through the window I\ at 
time ti, J 2 at time t 2 and so on. The windows I\, ■ ■ ■ ,I n -i are specified by the cylindrical 
interval J[N~] = J(N~ ) x x • ■ ■ x i x R T \ N . Since all the path in the 

integral (1111) have the right end fixed, there is no integration over the last instance t. We 
use the notation N~ = {t\,... ,t n ~ 1 } for the list of points in time excluding the last one. 

Note, that G is not probability function in the classical sense. Based on the rules of 
quantum mechanics, we state that it is reasonable to introduce “probability distribution 
functions” which take on complex values. Also we would like to remark, that the solution 
to the Schrodinger equation given by theorem [2] possesses a property of a Green function: 

lim ££[«(*,.-)] = S(x). 

Since the paths xt = x : [0,f] -ft have fixed both left and right ends, it follows from the 
definition of the marginal expectation E* t [U(xT-)\ that if t tends to zero, then E* t [U(xT-)\ 
is zero unless x = 0. 

Also, the solution E* t [U(xT-)\ is purely local, it is defined only in a space-time neigh¬ 
borhood of (x, t). In the next section we extend it uniquely to a global solution and give 
it representation in an operator form. First, let us introduce some notation: denote by 
K a compact set. Let T>k be a subset of C'^°(M)-functions that have their support in K : 
Dk — {<f> G C“(R)|supp <f C A'}; let V = U K V K . 

Corollary 4. Let K c M &e a compact set, 0 < t < t be fixed, f(x) G V K , V(y) be 
continuous at y = x G K, then there exist a neighborhood of (x,t) G K x (0, r] where the 
function 

-fi loc (x,t)= f E* t [U(x T -)]f(y)dy (13) 

J R 
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is a solution to the following IVP 

dif(x, t ) 


= i 


irk _ l VM 

2 dx 2 2 1 ' 


i>(x,t), 


dt 

*l>(x,Q) = f(x). 

Proof. First, we show that if loc (x,t) satisfies the initial condition (1151) : 

V ^ CM ) u = ( [ eiMvt-)] 


(14) 

(15) 


t =o 


f(y)dyj 


= lim 

t->o 


E*tlht(y T -)\f(y)dy 


As t tends to zero the condition on the integrand becomes x(t) —> x(0), thus the marginal 
expectation E* t yields one if and only if x=y and otherwise it yields zero, that means 

^ loc (x,t) I, „= [ E* t [U(y T -)]\ t=0 f(y)dy 


= / S(x-y)f(y)dy 

J R 

= f(x)- 

Thus the initial condition is satisfied. Next we prove that if loc (x,t) solves the equation (U4l) . 
We use the fact that if loc (x,t) is a locally integrable function. We have 




\v 


= / l 


i d 2 


2 dx 2 


KtM x T-)]f(y)dy 

V ) E* t [U(x T -)]f(y)dy 


d_ 

dt 


d_ 

dt 


E* xt \U(x T -)\f{y)dy by theorem [2] 
E*tP{x T -)\f(y)dy 


= 4p«( x , t) . 


□ 


4 GLOBAL SOLUTION TO THE SCHRODINGER EQUATION 
4.1 Class of potentials <SA(M) 

Our goal is to compare the solution for the Schrodinger equation via Henstock (fl3|) to the 
solution via Trotter’s theorem in the form of a one-parameter unitary group (J8J) . Since the 
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convergence of the limit in Trotter’s formula OH) requires the Hamiltonian to be essentially 
selfadjoint, we limit ourselves to considering a special class of functions serving as potentials. 


Definition 11. Let V : R —> R be a real measurable function. We say that V belongs to the 
class dh4(M) if and only if the operator (—A + V) is essentially selfadjoint on L 2 ' 


The question of giving precise mathematical conditions for essential selfadjointness of 
physical Schrodinger operators (—A + V) has a long history going back to the two pioneering 


papers of T. Kato 


2a 


30] in 1951. In a formulation given by J. von Neumann, the question 


is about the distinction between what in the math physics literature is called formally 
selfadjoint (alias, Hermitian, or symmetric), as opposed to selfadjoint. Typically the initial 
dense domain where the given Schrodinger operator (—A + V) is naturally defined needs to 
be completed in a graph closure, yielding thus the closure of (—A + V), as a closed Hermitian 
operator. If the closure is selfadjoint, we say that (—A + V) is essentially selfadjoint. In 
general settings, it may not be. But it is the latter selfadjointness condition which allows 
us to write a spectral resolution for (—A + V), and also to create effective path space 
representations for the associated Schrodinger equation. 

By J. von Neumann’s theory, we know that (—A + V) is essentially selfadjoint if and 
only if it has indices (0, 0). While the von Neumann indices have a physical interpretation, 
to decide essential selfadjointness is always a subtle question, both from the standpoint 
of mathematics and physics. We are interested in effective path space representations for 
the Schrodinger equation. But since we rely on first having a resolution to the essential 


selfadjointness question, we include here the following citations 


29 


35]. Each of these papers, 


starting with Kato 1951, offers useful insight into the question. For concrete examples we 
give the following conditions on potential function V that guarantee that it falls into the 
class dh4.(M): 


1. V(x) >0 and lim^i^oo V(x) = oo 

2 . —A + V has a dense set of analytic vectors [37]. 


3. For a system of N + 1 particles in 3-dimensional space denote the position vector of 
ith particle by r t = ( x^y^Zi ), then V can be expressed as 


N 


1 ,N 


v (n, • • • , 7-jv) = V'(r i, • • • , t N ) + Vijin - rj), 


1=1 


i<j 
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where V' is bounded in the whole configuration space, and Vij (0 < i < j < s), defined 
in 3-dimensional space, are locally square-integrable and bounded at infinity (see 
In other words, there exist two constants C and R such that 


2 *. 


' r<R 


\V'(r ir -- ,r N )\<C, 

\Vij(x,y, z)\ 2 dxdydz < C 2 , (0 < i < j < N) 

| Vij(x, y,z)\ < C (r > R,r = (x 2 + y 2 + z 2 ) 1 ^ 2 , 0 <i<j< N ). 


An example of the potential for which (—A + V ) is not essentially selfadjoint on L 2 
is V{x) = —x 4 (see [10]). The reason is that the classical motion reaches ±oo in finite time 


fix an energy E > 0 then (A) 2 — x 4 = E yields 


too 


dx 


yjE + x 4 


< oo. 


4.2 Operator JCt(x,f(-)) 

Let 0 G V, {Oi} be a finite open cover of supp (j) and /% be a partition of unity subordinate 
to {Oi}, i.e. the following conditions are satisfied: (1) f3i > 0, (2) j3i G C^°(M), (3) supp /3i C 
Oi, (4) = 1 for all x G supp Then we can write (ft = Pi<ft and define a linear 

map JCt(x, •) : T>k —> V as follows 

:= [ E* t [U(x T -)]f(y)dy 

Jr 

with an action on a test function «ft(x) G T>: 

< £t(x, /(■))> 0 (®) > = £t(x, f(-)),/3i<ft(x) > . (16) 

In the context of the corollary HJ K, t (x,f(-)) exists only locally in x, so to define its action 
as a distribution on D, we use the partition of unity. Representing <ft> = Yh A0 allows us to 
pair IC t (x, /(•)) with /3i<ft(x) that are non-zero in a sufficiently small neighborhood that can 
be chosen according to the corollary |4j Thus using linearity we define < JC t (x, /(•)), <ft>{x) > 
as a sum of pairing of K, t (x, /(•)) with j3i(ft{x), so the expression (TT61) makes sense. 

Next corollary extends the domain of fC t (x,f(-)) to L 2 (M). 
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Corollary 5. Under the conditions in the statement of the corollary with f(x) G L 2 (M) 
the function 


° c (x, t) = /C t (x,f(-)) 


is a local solution to the IVP p4\)-[T5\). 


(17) 


Proof. Let f(x) G L 2 (R), then we can write f n (x) 



and 


f(x) = lirn f n (x). 

n—too 


For all n G N, [—n,n] is compact, thus JC t (x, /„(•)) is well defined and by corollary []] is 
a local solution to (fl4l) - (fl5l) with the initial condition f>(x, 0) = / n (x). By interchanging 
limits with the integral sign in the proof of the corollary |4l we can show that the solution 
to (1T4|) - (TT5|) for f(x) G L 2 (M) is obtained as follows: 


ip loc (x,t ) = K. t (x,f(-)) = lirn lC t (x, /„(•))■ 

n—>oo 

□ 


4.3 Main result 


Theorem 6. If V(x) is a continuous function in then each of the local solutions to 

the Schrodinger equation f> loc (x,t) given by corollary^ has a unique extension if 9lob (x,t) for 
allt&R and iGl. The extension fj 9lob (x, t) solves the following IVP with f G L 2 (M): 


d'ip 9lob (x, t) 

dt 

iP 9lob (x , 0) 


. n d 2 

2 dx 2 



f(x), iGl 




tGl,a:Gl, 


(18) 

(19) 


Moreover, f/ 3loh {x, t) is represented by a unitary one-parameter strongly continuous group in 
L 2 (R): 

r'°\x,t) = K t (x,f(-)). (20) 


Proof. For V G 5^4(M) the operator —A + V is essentially selfadjoint, thus by Stone’s 


theorem 38] the solution to (1T8|) - (TT9|) is given as a strongly continuous group of unitary 
operators if(x,t) = {U t f){x) acting in L 2 (R). The operator U 1 is related to the operators 
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K f and UP fl6l)- (17|l via Trotter’s limit in the strong operator topology 


0 


V’OM) = 77)7) 


= lim (K^MfX f(x) 

n—>oo \ J 


= lim 


n->oo \2nit 


n \ 2 


n 

> 3 =1 


e [i" Xj- 1| 2 ** V(^)] f( Xn ^ Xl '"dx r 


n times 


= lim 


n—>oo V 27TZt 


U \ 2 




m I x j x J-ll~ 

2 (tT^p 


' f(x n )dxi... dx n . 


n times 


By corollary |5] there exist a local solution to the IVP (TT8|) - (TT9|) : 7 oc 7> t) = K t (x,f(-)), 
for small f. We would like to compare the local solution ijj loc (x,t) and the global Trotter 
solution ([/*/) (;r). 


Lemma 7. Under the conditions in the statement of the theorem for any f®<f>E L 2 (Jst)®T> 
there exist e > 0 such that for all |£| < e 


< (7/) 7)77) >=< Kt{x, f{-)),<t>{x) > . 


Proof. By the Schwartz kernel theorem (e.g. see theorem 5.2.1 in 39]) for the linear map 
/Q7, •) : L 2 (M) —> V defined in (fT6|) . there is one and only one distribution K on L 2 (M) ®V 
such that 


< 1Ct{x, /(■))> 07) >= K t(f ® 0) 


for all functions f(x) E L 2 (M), 07 ) E T>. Then the following is well-defined for |s| < |t| < e 


KsiiU^f) ® 0) =< K s {x, U f - S f), 0 > 

(< k , o 7 , f/ 7 (-)), 07 ) >=< (f/7) 7)> 07) >, if s = o 

1 < /Ct7,/(•)), 07) >, ifs = t. 
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Differentiate this expression in variable s, using the corollary [5j 

d ^ r ( TJ t-s t \ r < tC s+h (x,U t - s ~ h f),(j)>- 

— </C 8 (x,U /),</>> = hm- 7 

CIS h—> 0 iX 

<JC s+h (x , U t — h f),(j>>- 


= lim 

h—>0 


<K s {x ) U t ~ s f)A> 

<JC s (x,U t - s - h f),cl>> 


, < /C s (x, £/“/), </> > - < /C s (x, 0 > 

+ lim---- 

/i —h 

= i[< (A -V)K.(x, (U‘-‘f)),4>> 

+ (< — (A - V) tC s (x, (U‘-f)),4>>)] 

= 0 . 


Since the derivative is constant, we arrive to the equality for all f(x) G L 2 (M), (j>{x) G T> 
and \t\ < e\ 

< (U*f) (x),<l>(x) >=< JC t {xJ(-)),(j)(x) > . 

□ 


Since the operator [/*(•) coincides with K t (x, •) for small \t\ < £, where e is from lemma 
[3 we get K t (x, /(•)) G Lf oc (M). Also ^ 4 (') is unitary and possesses the semigroup property, 
thus K t (x, •) is also unitary and locally has the semigroup property 

K t (x, K s (x, •)) = K t+S (x, •) (21) 


for |i|, |s|, \t + s| < £. Now we can extend the local solution JC t (x ,/(•)) in time. For any 
i G I, there exist an integer m such that \t/m\ < e. Then by the semigroup property we 
have 


£t(x, /(•)) = /(-))) m = Kt /m {x, • • • lC t/m (x, /(•)) ■ ■ ■ )■ 

Local semigroup (|2B and unitary properties allow us to extend JC t (x,f(-)) globally as a 
one-parameter unitary group. To carry out the steps for going from local semigroup to a 


4Q 


44], JC t (x, /(•)) is also strongly continuous w.r.t. 


global unitary one-parameter group see 
L 2 since D 4 possesses this property. 

Hence, in conclusion, we get the global solution to the Schrodinger equation (IT8l) - (fT9l) 
realized as a Muldowney path-space average via its local ”complex integrals.” This solution 
is a strongly continuous group of unitary operators acting in L 2 ‘ 


r l °\x,t ) := JCt(x, f (•)). 


□ 
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5. SUMMARY 


As a summary, we give the outline of the technical steps in the proof of the main result. 
Because of the assumed essential selfadjointness, we automatically have the existence of the 
unitary one-parameter group U t from the spectral theorem, so a solution to the Schrodinger 
equation as a strongly continuous group of unitary operators U l acting in L 2 (M). And this 
unitary one-parameter group U l is therefore realized as a Trotter approximation. The latter 
allows us in turn to also realize it locally as Muldowney path-space average /Q. 

However, the drawback is that while the latter ]C t arises from the average operation 
based on Feynman’s intuitive idea, and on the rigorous Muldowney path-space calculus, 
it does not yield operators in L 2 (M). This is because it is based on path-space averages 
with complex densities. So /C t only yields a local solution to the Schrodinger equation as 
a system of Schwartz-distribution kernels, and so it does not yield operators in the Hilbert 
space L 2 (M). The L 2 realization exists only via the unitary one parameter gro up U t ) so the 
two must be connected in a precise way. This requires the machinery from 40-1441] which in 
turn allows a precise linking from local to global. 

So our piecing together the local path-space formulas to a global solution to the 
Schrodinger equation is accomplished in the following sequence of steps, and with us making 
use of the local-to-global theory from 


40 


44] as follows. We first show that: (i) U t and 


JC t agree locally on test functions, (ii) Using the one-parameter group-property for [A , we 
conclude that K t satisfies a local semigroup property, (iii) Since U t is unitary in L 2 for all 
t in M, we conclude that fC t extends canonically from initially being only a local semigroup, 
to a global realization as a unitary one-parameter group, (iv) Strong continuity w.r.t. L 2 is 
automatic since U 1 is already strongly continuous. 

In conclusion, we would like to remark that the results of the main theorem [b] hold in 
more general settings for any potential V for which the Trotter’s limit is convergent flHJ). 
The restriction V to the class iSM(M) guarantees the convergence of the Trotter’s limit, but 
it’s not an equivalent condition. 
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